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INTRODUCTION 


The  partially  plastic  problem  of  pressurized  thick-walled  cylinder  Is  of 
practical  Importance  to  pressure  vessels  and  the  autofrettage  process  of  gun 
barrels.  Many  solutions  for  this  problem  have  been  reported  (refs  1-7). 

For  thick  tubes  under  very  high  pressure  operation,  the  elastic-plastic 
material  model  should  be  represented  by  the  von  Mises'  yield  criterion, 
Prandtl-Reuss'  Incremental  stress-strain  laws,  the  strain-hardening,  and 
compressibility  (ref  8).  However,  a  closed-form  solution  exists  only  In  the 
plane-strain  case  neglecting  strain-hardening  and  compressibility. 

For  the  generalized  plane-strain  problems  considered  here,  numerical 
solutions  were  reported  by  the  finite-difference  method  (refs  4,7)  and  finite- 
element  method  (ref  5) .  The  incremental  displacements  were  used  as  the  basic 
unknowns  and  a  displacement  function  was  assisted  in  the  finite-element  method 
(ref  5).  The  Incremental  stresses  and  strains  were  used  in  Reference  4  as  the 
basic  unknowns,  but  only  the  Incremental  strains  were  used  in  Reference  7. 

The  spatial  discretization  used  in  References  4  and  7  was  based  on  the  forward 
difference  scheme  and  a  fixed  sequence  of  Incremental  loading  was  used. 

In  this  report,  a  new  method  is  developed  and  more  accurate  numerical 
results  are  obtained.  The  incremental  strains  are  chosen  as  the  basic 
unknowns  in  the  finite-difference  formulation.  Both  strain-hardening  and 
ideally-plastic  materials  can  be  considered.  The  spatial  discretization  is 
based  on  the  central  difference  scheme  and  the  Incremental  sizes  of  the 
applied  loading  are  determined  automatically  in  the  program.  The  incremental 


References  are  listed  at  the  end  of  this  report. 


results  are  calculated  directly  and  no  iteration  is  needed.  The  convergence 
of  the  approach  will  be  discussed  and  more  accurate  results  will  be  reported. 


BASIC  EQUATIONS 

Assuming  small  strain  and  no  body  forces  in  the  axlsymmetric  state  of 
generalized  plane-strain,  the  radial  and  tangential  stresses,  or  and  09,  must 
satisfy  the  equilibrium  equation, 

r(3cr/3r)  ■  09  -  or  (1) 

and  the  corresponding  strains,  and  eg,  are  given  in  terms  of  the  radial 
displacement,  u,  by 

er  -  3u/3r  ,  eg  ■  u/r  (2) 

It  follows  that  the  strains  must  satisfy  the  equation  of  compatibility 

r(3eg/3r)  -  er  -  eg  (3) 

If  the  material  is  assumed  to  be  elastic-plastic,  obeying  the  Mises*  yield 
criterion,  the  Prandtl-Reuss  flow  theory,  and  the  isotropic  hardening  law,  the 
stress-strain  relations  are  (ref  1): 


-  doi'/2G  +  (l/Doi'&o/ltiL') 

(M 

do  >  0  for  i  “  r, 0,z 

de*  -  E“  l(l-2  v)dom 

(5) 

where  E,  v  are  Young's  modulus,  Poisson's  ratio,  respectively, 

2G  -  E/(l+v), 

%  -  (Er+egfez)/3  ,  et'  -  e*  -  % 

°m  *  (or+09+oz)/3  ,  ai'  -  -  On 

a  -  (l//?)((ar-og)2  +  (cg-az) 2  +  (ffz-<*r)2]1/2  >  <*>  (6) 


and  0o  is  the  yield  stress  in  simple  tension  or  compression.  For  a  strain- 
hardening  material,  H'  is  the  slope  of  the  effective  stress/plastic  strain 


curve 


0  -  H(  JdeP) 


(7) 


For  an  ldeally-plastic.  material  (H'  ■  0),  the  quantity  (3/2)do/ (oH* )  is 
replaced  by  dX,  a  positive  factor  of  proportionality.  When  a  <  ff0  or  da  <  0, 
the  state  of  stress  is  elastic  and  the  second  term  in  Eq.  (4)  disappears. 
Following  Yamada  et  al  (ref  9),  Eqs.  (4)  and  (5)  can  be  rewritten  in  an 
incremental  form 

d«i  ■  d{4d&4  for  i,j  ■  r,9,z 

(8) 

and 

d^j  *  2G[v/(l-2v)  +  -  Oj'Oj'/S] 

where 

2  1  , 

S  •  -  (1  +  -  H'/G)o2  ,  H’/E  -  u/(l-u»)  (9) 

3  3 

uE  is  the  slope  of  the  effective  stress-strain  curve,  and  is  the  Kronecker 
delta. 

This  form  was  used  in  the  finite-element  formulation  for  solving  elastic- 
plastic  thick-walled  tube  problems  (ref  5).  In  the  following  section,  the 
Incremental  stress-strain  matrix  will  be  used  in  the  finite-difference 
formulation. 

FINITE-DIFFERENCE  FORMULATION 

Consider  a  thick-walled  cylinder  of  inner  radius  a  and  external  radius  b. 
The  tube  is  subjected  to  inner  pressure  p,  external  pressure  q,  and  end  force 
f.  The  elastic  solution  for  this  problem  is  well-known  and  the  pressure  p*, 
q*,  or  f*  required  to  cause  initial  yielding  can  be  determined  by  using  the 
Mises'  yield  criterion.  For  loading  beyond  the  elastic  limit,  an  incremental 
approach  of  the  finite-difference  formulation  is  used.  The  cross-section  of 

3 


the  tube  is.  divided  into  n  rings  with  r^-a,^ . r^-p . rn+l*b  where  p  is 

the  radius  of  the  elastic-plastic  interface.  At  the  beginning  of  each 
Increment  of  loading,  the  distribution  of  displacements,  strains,  and  stresses 
are  assumed  to  be  known  and  we  want  to  determine  Au,  ACj,  Aeg,  Acz,  Aor,  Aog, 
Aaz  at  all  grid  points.  Since  the  incremental  stresses  are  related  to  the 
Incremental  strains 'by  the  incremental  form  (Eq.  (8))  and  Au  *  rAeg,  there 
exists  only  three  unknowns  at  each  station  that  have  to  be  determined  for  each 
Increment  of  loading.  Accounting  for  the  fact  that  the  axial  strain  ez  is 
independent  of  r,  the  unknown  variables  in  the  present  formulation  are  (Aeg)^, 
(Aer)^,  for  i  ■  1,2, . . .  ,n,n+l,  and  Aez. 

The  equation  of  equilibrium  (1)  and  the  equation  of  compatibility  (3)  are 
valid  for  both  the  elastic  and  the  plastic  regions  of  a  thick-walled  tube.  A 
first-order-correct  finite-difference  analog  of  these  two  equations  at  i  * 

1 . .  has  been  given  in  References  4  and  7.  Other  finite-difference  forms 

can  be  written.  In  this  report,  the  difference  equations  given  below  are 
second-order-correct.  The  equation  of  compatibility  (3)  and  equation  of 
equilibrium  (1)  are  replaced,  respectively,  by 


and 

where 


cli( Aeg)j[  +  c2i(Aer)i  +  C3i(Aeg)i+i  +  c4ji(Aer)i+i  -  cs* 

cli(Aor)i  +  C2i(Aog)i  +  C4i(Aor)i+i  +  C4i(Aog)i+i  -  eg* 


(10) 

(ID 


c5i  “  Cli(e6>i  -  C2i(er)i  -  C3i(eg)i+1  -  C4i(er)i+i 
C61  “  "  cii(or)1  -  C2i(ae)i  “  c3i( °r)i+l  “  C4i(og)i+i 


and 

?i  -  ri+i/ri  (12) 

With  the  aid  of  the  Incremental  stress-strain  relations  (8),  Eq.  (11)  can 
be  written  as 

c7i(Ae8)i  +  c8i(Aer)i  +  C9t(deQ)i+i  +  cioi(Acr)i+i  +  ciu(Aez)  -  c^i  (13) 
where 

c7i  ■  cu(di2)i  +  C2i(d22)i  »  C8t  -  cii(dn)i  +  C2i(d2l)i 

C91  "  C3i(dl2)i+1  +  C4i(d22)l+1  .  CIO!  “  C31<<ill)l+1  +  C4i(d22>i+1 

clli  ■  cu(di3)i  +  C2i(d23)l  +  C3i(d13)i+1  +  C4i(d23)i+1  (14) 

The  boundary  conditions  for  the  problem  are 

Aor(a»t)  -  -Ap  ,  Aor(b,t)  -  -Aq 
n 

ir  l  [r^Ao^  +  ri+i(Aoz)1+i](ri+i-ri)  -  jnra2Ap  +  Af  (15) 

1-1 

where  y  Is  zero  for  open-end  tubes,  and  one  for  closed-end  tubes.  Using  the 
Incremental  relations  (8),  we  rewrite  Eq.  (15)  as 

(di2)i(Aea)i  +  (dj.i)l  (ACy)!  +  (di3)iAez  «  -Ap  (Id) 

<<*12)n+l(Aee)iri-l  +  (dll)n+-l(Aer)n+l  +  (d13>n+lAez  *  (17) 

and 

n 

l  lci2l(Aee)i  +  c13i( Aer)i  +  ci4i(A6e)i+i  +  ci5i(AEr)1+1 
1-1 

+  c16i(  Aez)  1  *  Ua2Ap  +  Af/  if  (18) 

where 

c12i  *  (ri+l-ri)rt(d32)i  ,  ci3i  ■  Ui+l“ri)ri(d3l)i 
c14i  *  (ri+l"ri)ri+l(d32)i+l  ,  cisi  "  (ri+l“*i)ri+l(d3l)l+l 

c16i  -  (ri+l“ri) tri(d33)1  +  ri+i(d33)1+l]  (19) 

Now  we  can  form  a  system  of  2n+3  equations  for  solving  2n+3  unknowns,  (AeQ)^, 
(Aer)i,  at  1  -  1,2, ... ,n,tt+l  and  Aez.  Equations  (16),  (17),  and  (18)  are 


taken  as  the  first  and  last  two  equations,  respectively,  and  the  other  2n 
equations  are  set  up  at  i  ■  l,2,...,n  using  Eqs.  (10)  and  (13).  The  final 
system  is  an  unsymmetric  matrix  of  arrow  type  with  the  nonzero  terms  appearing 
in  the  last  row  and  column  and  others  clustering  about  the  main  diagonal,  two 
below  and  two  above.  In  the  computer  program  which  was  developed,  the 
dimensionless  quantities  r/a,  Zex/a0,  Eeg/Og,  Eez/o0,  ar/o0,  ag/aQ,  oz/o0, 
p/°o*  ^ao»  ^/(ir®2°o)  were  used  in  the  formulation  and  the  Gaussin  elimination 
method  was  used  to  solve  these  equations.  All  calculations  were  carried  out 
on  IBM  4341  with  double  precision  to  reduce  round-off  errors. 

OPTIMAL  INCREMENTAL  LOADING 

Given  any  combination  of  incremental-loading  (Ap,  Aq,  or  Af ) ,  we  can  now 
determine  all  Incremental  results  (displacements,  strains,  and  stresses) 
directly.  No  iteration  Is  needed,  while  in  Reference  4,  many  iterations  in 
each  step  were  required  because  a  value  for  Aez  was  assumed.  The  sizes  of 
incremental- loading  should  be  chosen  properly  in  order  to  obtain  accurate 
results  at  a  reasonable  cost.  When  the  total  applied  pressure  p  is  given,  it 
is  natural  to  divide  the  loading  path  in  m  equal  fixed  increments  such  as  Ap  - 
(p-p*)/a.  Larger  values  of  m  give  more  accurate  results.  A  sequence  of 
decreasing  load-increments  is  a  better  choice  than  that  of  equal  increments. 

In  order  to  increase  the  efficiency  without  affecting  the  accuracy,  an 
adaptive  algorithm  has  been  implemented  on  the  basis  of  a  scaled  incremental- 
loading  approach  (ref  5). 

In  each  step,  a  dummy  load-increment  such  as  Ap  is  applied  and  the 
incremental  results  Ao^  for  i  -  r,0,z  at  all  grids  are  determined.  For  all 


grid  points  at  which  o  ■  llo^M  <  o,  we  compute  the  scaler  o’ s  by  the  formula 
a  *  -  (T  +  [  T2  +  4|  UaJ  |  2<0o*-|  I  Oil  I  2)] 1/2  }/ 1  UaJ  |2  (20) 

where 

r  -  I  I  a± 1 1  2  +  1  |  A<?±  |  !  2  -  1 1  Ci+dOi  |  |  2  (21) 

and  ||ajj|,  | | Ao£ | | ,  ||a^+  Aa^ | |  are  computed  by 

I  |<Jil  |2  -  -  [(or-<J0)2  +  (ag-az)2  +  (oz-or)2]  (22) 

Let  X  be  the  minimum  of  the  a’ s.  Then  X  is  the  load-increment  factor  just 
sufficient  to  yield  one  additional  point.  A  sequence  of  X(j)  can  be 
determined  for  all  steps  j  -  l,2,...m  and  the  updated  results  are 

p(j)  -  p(j-l)  +  xC ApC J) 

+  xO)do^Cj)  ,  etc.  (23) 

This  sequence  of  Incremental  loading  is  optimal  for  the  present  problem 
because  all  the  coefficients  c's  in  Eqs.  (12),  (14),  and  (19)  are  functions  of 
the  previous  stresses  and  strains  •• 

CONVERGENCE  STUDY 

In  order  to  demonstrate  the  accuracy  of  the  approach,  four  convergence 
studies  are  made.  Consider  a  thick-walled  tube  of  wall  ratio  b/a  ■  2  and 
subjected  to  internal  pressure  only.  The  cross-section  of  the  tube  is  divided 
into  n  rings  of  equal  thicknesses,  i.e.,  h  ■  (b-a)/n.  The  first  problem  is  a 
closed-end  tube  loaded  in  the  elastic  range  with  G  »  10s/3  psi,  v  -  0.3,  p  ■  5 
psi.  The  numerical  results  with  n  -  10,  20,  50,  100  are  shown  in  Table  I 
together  with  the  Lame"  solution  for  the  hoop  stresses  and  strains  at  the 
boundaries  a  and  b.  The  numerical  results  are  correct  up  to  four  digits  with 
n  -  100. 
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TABLE  X.  ELASTIC  SOLUTION  FOR  A  CLOSED-END  TUBE 
(b/a  »  2,  G  -  105/3  psi,  v  -  0.3,  p  ■  5  psi) 


°0|a  !  a9|b  !  Ee6|a  i  Ee0|b  | 


I 

8.3333  |  3.3333  |  9.3333  I  2.8333  I 


10 

1 

1 

8.3000 

1 

1 

3.3000  | 

3.3250  | 

20 

1 

1 

I 

8.3256 

1 

1 

1 

50 

1 

I 

1 

8.3320 

1 

1 

1 

3.3320  | 

100 

1 

1 

8.3330 

1 

1 

3.3330  | 

9.3000 

9.3250 

9.3320 

9.3330 


I  2.8333  I 


I  2.8000  | 

I  I 

I  2.8250  | 

I  I 

I  2.8320  ! 

I  I 

I  2.8330  | 


The  second  problem  Is  the  Initial  yielding  solution  for  a  plane-strain 
tube  with  E/o0  ■  200,  v  *  0.3,  ez  *  0.  The  numerical  results  with  n  ■  10,  20, 
50,  100,  200,  are  shown  in  Table  II  together  with  the  exact  solution  for  the 
dimensionless  p  -  p/<J0,  o0  -  a0/oo,  e9  *  (E/<Jo)(e0)  at  r  »  a  and  b. 

TABLE  II.  INITIAL  YIELDING  SOLUTION  FOR  A  PLANE-STRAIN  TUBE 
(b/a  -  2,  E/o0  -  200,  v  -  0.3) 

]  i  i  :  i  :  r~;  i  : - 1  :  r 


Exact 


P 

°9|a  ] 

°9|b 

e9|a 

43229 

.72049  | 

.28820 

.82424 

00110 

-.00107  | 

-.00216 

-.00054 

00028 

-.00027  j 

-.00054 

-.00014 

00005 

-.00004  | 

-.00009 

-.00002 

00001 

-a  00001  | 

j 

-.00003 

-.00001 

00001 

.00000  1 
- _ 1, 

-.00001 

.00000 

( — r  r 

I  e9|b  I 


.2622 

-.00065 

-.00016 

-.00003 

-.00001 

.00000 


The  third  problem  is  the  elastic-perfectly  plastic  solution  for  a  plane- 


strain  tube  with  b/a  -  2,  E/c0  -  200,  v  -  0.3,  u  -  0,  ez  -  0.  The  numerical 
results  with  n  -  10,  20,  50,  100,  200  are  shown  in  Table  III  for  the  pressure 
and  displacement  at  the  bore  corresponding  to  50  percent  and  100  percent 
overstrain,  i.e.,  p/a  ■  1.5  and  2.0. 

TABLE  III.  ELASTIC-PERFECTLY  PLASTIC  SOLUTION  FOR  A  PLANE-STRAIN  TUBE 


(b/a  -  2,  E/00  -  200,  v  -  0.3) 


1 - 

1 

1 

1 

1 

1 

50% 

Overstrain 

~T 

1 

! 

T 

100%  Overstrain  1 

1 

1 

I  n 

1 

1 

1 

1 

P 

r 

i 

i 

e0l  a 

1 

1 

P 

|  e0la  j 

1 

1  10 

1 

1 

1 

1 

.71863 

i 

i 

i 

1.97725 

1 

1 

i 

.79849 

I  3.68937  | 

1 

1 

I  20 

1 

1 

1 

| 

.71825 

1 

i 

1 

1.96969 

1 

1 

1 

.79790 

1  1 

1  3.66831  1 

1  1 

1  50 

1 

1 

1 

.71808 

I 

i 

i 

1.96780 

1 

1 

i 

.79760 

I  1 

1  3.66280  | 

1 

I  100 

1 

1 

1 

.71803 

i 

1.96759 

1 

i 

.79751 

1  1 

1  3.66211  | 

t  1 

I  200 

1 _ 

1 

1 

.71801 

1 

i 

1.96757 

1 

1 

.79747 

i  1 

1  3.66199  I 

1  1 

The  numerical  results  converge  and  are  accurate  up  to  four  digits  with  n  ■ 
100.  There  is  no  closed-form  solution  available  for  comparison.  The  famous 
paper  by  Hodge  and  White  (ref  2)  has  been  used  quite  often  as  a  basis  for 
assessing  the  accuracy  of  other  approximate  methods.  However,  the  numerical 
Integration  is  cumbersome.  The  present  formulation  is  much  simpler  and  seems 
more  accurate.  In  order  to  further  demonstrate  the  accuracy  of  the  present 
approach,  a  convergence  study  for  an  incompressible,  ideally-plastic  thick 
tube  in  plane-strain  condition  has  been  made  and  compared  with  exact  solution 
(ref  2).  The  numerical  results  for  a  nearly  incompressible  material  (v  ■ 


0.49999)  are  shown  in  Table  IV  together  with  the  exact  solution  (v  ■  1/2)  for 
the  internal  pressure  and  the  displacement  at  the  bore  corresponding  to  50 
percent  and  100. percent  overstrain. 

TABLE  IV.  INCOMPRESSIBLE,  IDEALLY-PLASTIC  SOLUTION  FOR  A  PLANE-STRAIN  TUBE 


(b/a  « 

2,  E/o0  -  200, 

v  -  0.49999) 

1 - 

1 

,  . 

1 

1 

50%  Overstrain 

100* 

Overstrain  ! 

1 

1  . 

1 

1  n 

1 

P 

i  e0l  a 

P 

1 

1 

1 

e9la  1 

1 - 

1  10 

1 

.72069 

i  1.95714 

.80015 

r 

i 

i 

3.48806  I 

1 

r 

1 

I  20 

1 

.72077 

1  1.95072 

.80034 

l 

i 

i 

1 

3.47012  | 

1  50 

.72078 

S  1.94891 

.80038 

1 

i 

i 

3.46509  | 

1 

I  100 

1 

.72078 

I  1.94865 

.80038 

l 

i 

i 

1 

3.46437  | 

1 

1  Exact 

.72078 

I  1.94856 

.80038 

l 

i 

3.46410  | 

L 


We  may  thus  conclude  that  exact  solutions  can  be  obtained  by  this  numerical 
approach. 

ADDITIONAL  RESULTS 

After  establishing  the  convergence  and  accuracy  of  this  new  approach,  the 
numerical  results  for  more  general  problems  have  been  obtained.  Some  of  the 
additional  results  are  documented  here  for  future  comparison  by  others.  All 
numerical  results  presented  here  are  for  a  thick-walled  tube  with  wall  ratio 
b/a  ■  2,  E/o„  -  200,  v  *  0.3,  n  ■  100,  u>  *  Et/E  ■  0.1.  The  numerical  results 


are  accurate  up  to  four  or  five  digits.  Table  V  shows  the  results  of  the 
dimensionless  quantities  p/a0,  oe/o0,  °z/®o»  (E/Oq)***  (E/o0)ee,  (E/o0)ez  at 


2.0  in  a 


the  inside  or  elastic-plastic  boundary  p  for  p/a  *  1.0,  1.1,  1.2,..., 
plane-strain  tube  with  strain-hardening  parameter  <u  •  0.1.  Tables  VI  and  VII 
show  the  similar  results  of  the  dimensionless  stresses,  strains,  displacement 
at  the  bore  or  elastic-plastic  boundary  for  various  stages  of  elastic-plastic 
loadings  in  an  open-end  or  closed-end  tube,  respectively. 


TABLE  V.  ELASTIC-PLASTIC  SOLUTION  FOR  A  PLANE-STRAIN  TUBE 


(b/a  - 

2,  E/oq  -  200,  v  - 

0.3,  Et/E 

■ 

0.1,  n  ■ 

100,  ez  - 

0) 

p/a 

T 

1 

1 

1 

p/°o 

1  1 

I  Inside  1 

j  O0/Oo  j  O0/Oo|p 

T 

1  Inside 
j  °z/°o 

T 

1 

1 

1 

^  p 

r  E 

7  er|  a 

1  °o 

1 

1 

1 

1 

E  Ua 

°o  » 

1.0 

T 

1 

1 

.43229 

1 

|  .72049 

.72049 

I  .08646 

1 

T 

1 

1 

.08646 

I  -.67438 

T 

1 

| 

.82424 

1.1 

1 

1 

| 

.51296 

I  .66972 

l 

.75016 

1 

1  .05586 

j 

1 

1 

1 

.10453 

I  -.91636 

1 

1 

I 

1.00141 

1.2 

( 

1 

1 

.58362 

1 

I  .63049 

1 

.78249 

1 

I  .02537 

1 

1 

1 

I 

.12427 

|  -1.17381 

1 

1 

| 

1.20371 

1.3 

1 

1 

1 

.64556 

1 

I  .60225 
j 

.81739 

1 

I  -.00409 
j 

1 

1 

1 

.14566 

1  -1.44917 

1 

1 

1 

1.43022 

1.4 

1 

1 

i 

.69982 

I  .58427 

.85479 

1  -.03148 
| 

1 

1 

1 

.16866 

I  -1.73882 

1 

1 

f 

1.68000 

1.5 

1 

i 

.74725 

1  .57572 

.89457 

1  -.05589 

i 

1 

1 

1 

.19323 

I  -2.04304 

1 

1 

| 

1.95207 

1.6 

1 

1 

i 

.78851 

1  .57579 

.93667 

1 

1  -.07658 

1 

1 

1 

1 

.21932 

1  -2.36107 

1 

1 

1 

2.24535 

1.7 

1 

i 

.82416 

1  .58376 

1 

.98092 

1 

1  -.09302 

1 

1 

.24687 

1  -2.69205 

1 

1 

| 

2.55869 

1.8 

1 

1 

I 

.85466 

1 

1  .59898 

1.02718 

1 

I  -.10490 

1 

1 

1 

I 

.27581 

I  -3.03504 

1 

1 

1 

2.89081 

1.9 

1 

1 

i 

.88041 

1  .62086 

1.07530 

1 

1  -.11205 

1 

1 

1 

1 

.30606 

1  -3.38898 

1 

1 

1 

3.24034 

2.0 

1 

_L 

.90177 

1  .64885 

J _ 

1.12509 

1 

1  -.11447 

J _ 

1 

1 

_L 

.33753 

I  -3.75274 

J _ 

1 

1 

_L 

3.60578 

~T 


ELASTIC-PLASTIC  SOLUTION  FOR  A  CLOSED-END  TUBE 
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